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One of the shortcomings of the original theory of the geomorphologic unit hydrograph (GUH) is that 
it assumes that runoff is generated uniformly from the entire catchment area. It is now recognized that 
in many catchments much of the runoff during storm events is produced on partial areas which usually 
form on narrow bands along the stream network. A storm response model that includes runoff 
generation on partial areas by both Hortonian and Dunne mechanisms was recently developed by the 
authors. In this paper a methodology for integrating this partial area runoff generation model with the 
GUH-based runoff routing model is presented; this leads to a generalized GUH. The generalized GUH 
and the storm response model are then used to estimate physically based flood frequency distributions. 
In most previous work the initial moisture state of the catchment had been assumed to be constant for 
all the storms In this paper we relax this assumption and allow the initial moisture conditions to vary 
between storms. The resulting flood frequency distributions are cast in a scaled dimensionless 
framework where issues such as catchment scale and similarity can be conveniently addressed. A 
number of experiments are performed to study the sensitivity of the flood frequency response to some 
of the “similarity” parameters identified in this formulation. The results indicate that one of the most 
important components of the derived flood frequency model relates to the specification of processes 
within the runoff generation model; specifically the inclusion of both saturation excess (Dunne) and 
Horton infiltration excess runoff production mechanisms. The dominance of these mechanisms over- 
different return periods of the flood frequency distribution can significantly affect the distributional 
shape and confidence limits about the distribution. Comparisons with observed flood distributions 
seem to indicate that such mixed runoff production mechanisms influence flood distribution shape. 1 he 
sensitivity analysis also indicated that the incorporation of basin and rainfall storm scale also greatly 
influences the distributional shape of the flood frequency curve. 


1. Introduction 

Considerable research effort has been spent in recent 
years on the estimation of flood frequencies using the 
derived distribution approach. The first example of such an 
effort was described by Eagle son [1972]. Two more recent 
examples were presented by Hebson and Wood [1982] and 
Diaz-Granados el al. [1984], The rainfall-runofT models used 
by the latter are based on versions of the geomorphologic 
unit hydrograph (GUH) proposed by Rodriguez-Iturbe and 
Valdes [19791- 

Despite the simplicity and ease of parameterization of the 
GUH, a number of significant shortcomings still exist in the 
use of GUH-based models for the estimation of derived flood 
frequency. Recently, Moughamian et al. [1987] compared 
the approaches of Hebson and Wood [1982] and Diaz- 
Granados et al. [1984]. They found that both models per- 
formed poorly in every catchment studied when compared to 
sample distributions, suggesting that fundamental improve- 
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ments are needed before they can be applied with any 
confidence. 

The GUH of Rodriguez-Iturbe and Valdes [1979] is based 
on the assumption that rainfall excess is generated uniformly 
throughout the catchment area. In their application of the 
GUH to derived flood frequency estimation, both Hebson 
and Wood [1982] and Diaz-Granados et al. [1984] assumed 
two simple models of essentially Hortonian runoff genera- 
tion to calculate the rainfall excess. In reality, runoff gener- 
ation on catchments is much more complex. It is now 
recognized [Dunne, 1978] that in many catchments much of 
the runoff during storm events is produced on variable 
contributing areas which form narrow bands adjacent to the 
streams. The failure to incorporate these observed features 
of the catchment response into the GUH-based models is 
likely to be a source of major errors in the derived flood 
frequency estimates. 

In all of the examples cited above, the initial moisture 
condition of the catchment prior to the storm is assumed to 
be the same for all the storms. This is contrary to the reality 
in many catchments. Wood [1976] extended the flood fre- 
quency model of Eagleson [ 1972] in order to study the effects 
of parameter uncertainties. He found that uncertainty in the 
parameter that represented the initial moisture state of the 
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catchment could have a substantial effect on the predicted 
return periods. 

In addition, both Hebson and Wood [1982] and Diaz- 
Granados et al. [1984] used a lumped representation of the 
catchment response with spatially homogeneous rainfall 
based, however, on point rainfall statistics. This inconsis- 
tency gives rise to scale-dependent biases in the derived 
flood frequency distributions. Wood and Hebson [1986] 
overcame this problem by deriving a scale-independent flood 
frequency curve using an areal rainfall input distribution 
based upon areal rainfall similarity. 

Finally, Surkan [1969] has shown that catchments that are 
topologically similar (based for example on Horton’s order 
ratios) may yet produce different impulse response func- 
tions. This is the result of topological randomness whereby 
networks with different link-node configurations give rise to 
the same order ratios. For this reason, Beven [1986] has 
argued against generalizing the network for making predic- 
tions about specific catchments; instead Beven used a simple 
routing procedure based on a constant channel wave veloc- 
ity and network link histogram. More recently, Gupta et al. 
[1986] have presented an approach that utilizes the actual 
network structure directly for the routing of runoff through 
use of the so-called “width function.” This is a subject of 
intense current research interest. 

1.1. Scope of the Paper 

This is the third in a sequence of papers whose aim has 
been to provide a greater understanding of the interrelation- 
ships that underlie the storm response of catchments of 
different scales and physical characteristics by focusing on 
concepts of similarity. In the first paper, Wood and Hebson 
[1986] developed similarity relationships for flood frequency 
distributions that are independent of basin scale. They 
obtained a dimensionless flood frequency curve using the 
GUH basin response model of Rodriguez-Iturbe and Valdes 
[1979] under the assumptions of spatially homogeneous 
rainfalls and a simple Hortonian runoff generation model 
consisting of a constant contributing area during the storm 
events. 

Later, Sivapalan et al. [1987] relaxed these assumptions 
and developed a model of storm runoff generation due to 
spatially variable rainfalls on heterogeneous catchments 
taking account of the effects of catchment topography on the 
within-storm dynamics of runoff contributing areas. This 
model was also expressed in dimensionless form leading to 
the identification of five dimensionless catchment similarity 
parameters and three dimensionless auxiliary variables 
which govern the scaled storm response. 

This paper represents the logical next step toward the 
development of a physically based flood frequency curve. 
First, we generalize the existing GUH theory in order to 
incorporate, in a simple and parsimonious way, runoff gen- 
eration on partial areas by both the infiltration excess 
(Hortonian) and saturation excess (Dunne) mechanisms. The 
generalized GUH is structured in such a way that it can be 
coupled to the runoff generation model of Sivapalan et al. 
[1987] by means of the topography-soil index which is used 
to predict the areal distribution of soil moisture deficits and 
the proportion of contributing areas. This coupling permits 
the lumped runoff generation model of Sivapalan et al. 
[1987] to be integrated with the generalized GUH in a 
physically consistent manner. 


Second, the runoff generation model of Sivapalan et al. 
[1987] and the generalized GUH are combined together with 
assumed distributions of rainfall intensity and duration and 
an assumed distribution of the initial conditions to obtain the 
flood frequency curve. The scaled dimensionless formula- 
tions of Wood and Hebson [1986] and of Sivapalan et al. 
[1987] are essentially preserved in the derivations of the 
catchment storm response. The scaled flood frequency dis- 
tributions are used to address such issues as catchment scale 
effects and similarity. The sensitivities of the flood frequency 
curve to changes of various similarity parameters are also 
investigated. 

1.2. Outline of the Paper 

The paper begins with a summary of the GUH theory. The 
outline of the theory, given in section 2.1, is made suffi- 
ciently general so as to be able to incorporate partial area 
runoff generation. This is followed in section 2.2 by the 
introduction of a runoff generation model based on the 
topography-soil index. A methodology for combining the 
partial area runoff generation model with the GUH-based 
runoff routing model is presented next in sections 2.3 and 
2.4; this results in the generalized GUH. The sensitivities of 
the generalized GUH to two parameters of the partial area 
generation model are then studied. The results of this study 
are presented in section 2.5. 

Section 3 involves the estimation of the dimensionless 
flood frequency distribution. In section 3.1 the dimension- 
less peak discharge is calculated based on the dimensionless 
S curve derived from the scaled generalized GUH. Sections 
3.2 and 3.3 are devoted to the specification of frequency 
distributions of the scaled rainfall intensities, storm dura- 
tions, and the initial conditions. Section 3.4 describes the 
procedure for the estimation of the flood frequency distribu- 
tion and return period. Section 4 presents results of the 
sensitivity analyses carried out on the model with respect to 
a number of similarity parameters. 

2. Derivation of a Generalized GUH Based on 
Partial Area Runoff Generation 

2.1. A Review of the GUH 

Let T b denote the time of travel of any particle of surface 
runoff from the location of its generation to the catchment 
outlet. It is assumed that the time-variant instantaneous unit 
hydrograph (iuh) at the time r, /i(flr), is given by 

d 

H(t\T) = T P(T B ^t\T) (1) 

at 

where P( ) denotes the probability of the event given in the 
parentheses. Note that for all particles that are generated at 
time r, both T B and t are measured from the instant t. The 
right-hand side of (1) is simply the probability density 
function (pdf) of T B . 

Let fl be the highest order of the catchment stream 
network; C/(l < / < D) denotes a channel state of order i; and 
r,(l < i < fl) denotes an overland flow region or hillslope 
state of order We assume that runoff is generated only on 
the hillslopes; rainfall falling directly into stream channels is 
neglected. We can then define a collection S - {j} of paths s 
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which a particle of water may follow from the state r, where 
it is generated to the catchment outlet. Thus we have 

P(T b ^I\t)= X P(T s £t)P(s\r) (2) 

sSS 

where T s (s E S) is the travel time along a path 5 and P(s\t) 
is the probability of the given path s out of all paths 5, for 
given r. Neglecting the travel time over the hillslopes, any 
path s E S necessarily takes a form s = (x x , x 2 , m * * * x k >, 
where jc it * 2 , • * ' , x k E {c,-, i = 1, fl}. For a path s E S that 
originates in a hillslope of order / (note that x\ - C/), F(jIt) 
is given by 

P(s\r) = n,{T)P XlJt P IiXj -‘P Xl ^ t (3) 

where I1 /(t) is the probability that the particle of runoff came 
from a hillslope of order / and P xx is the transition proba- 
bility for the particle between order*, and Xj. In the general 
formulation it is equated to the ratio of the number of 
particles of runoff generated in all the hillslopes of order / to 
the total number of particles of runoff generated from the 
entire catchment. The generalized GUH is then obtained by 
combining (1), (2), and (3). The probabilities Ityr), for « = 1, 
in (3) are dependent on the runoff generation model. Their 
derivation is presented in section 2.3. The assumption ne- 
glecting hillslope travel times is reasonable for large catch- 
ments or narrow partial areas. The interaction between 
hillslope and channel travel times over a range of catchment 
scales is currently being investigated by the authors. 


2.2. Model of Runoff Generation 

The runoff generation model used in this paper is a 
simplified version of the conceptual model of Sivapaian et 
al. [1987]. It is based on the fundamental assumption that, 
under quasi-steady conditions, the difference between the 
local prestorm water table depth z x at a location * and its 
catchment-wide average z is linearly related to the corre- 
sponding difference between In (aT f /T x tan /3), a topography- 
soil index, and its catchment average A. This relationship can 
be expressed as 



In (4), a denotes the area draining through location * per unit 
contour length, T x is a local transmissivity parameter, and / 
is a hydrogeological constant for the catchment and is a 
measure of the decline of the saturated hydraulic conductiv- 
ity with depth. For soils for which this decline is exponen- 
tial, Seven [1986] showed that T x = K 0 /f where K Q is the 
saturated conductivity of the surface soil layer; A and T e are 
given by 



In T e = — f In T x dA 

A Ja 

Given z and the spatial pattern of values of the topogra- 
phy-soil index In (aT e IT x tan 0), equation (4) enables the 


prediction of the pattern of local initial water table depths for 
all points in the catchment. 

Generation of saturation excess runoff depends on the 
initial moisture storage deficit S x at any location x. For 
simplicity, we neglect the downslope redistribution of mois- 
ture within the duration of the storm and assume that 
saturation excess runoff will be generated wherever the 
cumulative infiltration M g {r) exceeds S x . It is assumed 
further that just before the storm the unsaturated zone 
moisture profile comes close to the case of complete gravity 
drainage; i.e., reaches “field capacity.” This assumption is 
consistent with field observations. S x may then be uniquely 
predicted from the depth to the water table. This unique 
functional relationship between z x and S x is denoted here by 
z x = y(5*). 

The initial contributing area A f (0) is obtained by using (4) 
to determine the value of the topography-soil index for 
which z x ^ ^ , where is the thickness of the capillary 
fringe and is assumed to be equal to the air-entry value of the 
matrix head for the soil. The contributing area is thus given 
by 

In ( — \ ^ A +/z -Jtyc 

\T X tan p J 

The contributing area A c expands with time during the 
storm, and at any time rafter the beginning of the storm, the 
dynamic contributing area is obtained by determining the 
locations where M g (r) ^ S x . In terms of the topography- soil 
index this condition can be expressed as 

In ( - aTt - ) * A + fl - fy[M g (r)] (6) 

\T X tan p) 

The model presented by Sivapaian et al. [1987] assumes 
that infiltration is controlled by the initial moisture content 
and saturated hydraulic conductivity of the surface layer and 
neglects the variation of these quantities with depth. This 
allows analytical solutions for the infiltration process to be 
developed. Both the rainfall intensity p and the surface 
hydraulic conductivity Kq are assumed to be spatially vari- 
able. Rainfall p is assumed to be gamma distributed within 
the catchment with mean p and coefficient of variation 
C vp . Kq is assumed to be lognormaily distributed with mean 
Rq and coefficient of variation C vK . The model uses quasi- 
analytical expressions for mean infiltration rate m g (T) and 
cumulative infiltration volume M g (r). These were derived by 
Sivapaian [1986] based on the Philip [1957] infiltration equa- 
tion and the time compression approximation. 

In this paper, a further simplification of the modeling of 
infiltration excess runoff is effected by approximating the 
infiltration process by a single lumped equation that uses the 
average moisture content over the noncontributing areas. 
This approximation is justified by the results of Sivapaian et 
al. [1987] which show that the variation of surface moisture 
content in space accounts for a negligible component of total 
runoff. 

Following Sivapaian et al. [1987], the model equations are 
expressed in dimensionless form in terms of five dimension- 
less similarity parameters and three auxiliary variables. The 
listing and definition of these eight parameters are presented 
in Appendix A. For more details, the reader is referred to 
Sivapaian et al. [1987]. 



46 


SlVAPALAN ET AL-: HYDROLOGIC SIMILARITY, 3 


2.3. Derivation of W^t) 

Let A be the area of the entire catchment. We denote by 
A* the total area of all the hillslopes that drain directly to 
streams of order i; 0, denotes the proportion of A* to A, i.e., 
0/ - A*/A and 10 f = 1. Expressions for 0/ for a third-order 
catchment in terms of Horton’s order ratios R A and R B are 
presented in Appendix B. 

Let A c (r) be the contributing area for the entire catchment 
defined as the area contributing direct runoff by the satura- 
tion excess mechanism. We also define by A*( t) the total 
contributing area from all the hillslopes which drain directly 
to streams of order /. Then we can easily establish the 
following relationship: 


A c (t) 

A 



(7) 


Let m f (r) be the areal mean infiltration rate over the 
non-contributing area of the catchment, A - A r (r). Here 
m g (r) is modeled using an average surface moisture content 
of the soil over the nonponded area. The areal average 
rainfall intensity, assumed constant in time during the storm 
event, is denoted by p. Then the rate of runoff generation at 
any time r from the entire catchment is given by 


, x Uc(t) 

<!(*) ~ A j — — p + 



[p - 


} 


( 8 ) 


We also make the scale approximation that the averages of 
rainfall and the infiltration rate taken over, say, all of the 
hillslopes that drain directly to streams of order i are equal to 
the catchment-wide means p and m g (r ), respectively. The 
mean runoff production rate from all the hillslopes that drain 
to streams in order i is then given by 


qfa) 


U»(r) 

-M-7F* + 



(9) 


The proportion II/(t) of the number of particles of runoff 
generated on all the hillslopes that drain to streams of order 
i f to the total number from the entire catchment is simply 
qfq\ thus we have 



At this point we approximate A*/A* in the infiltration 
excess component of (10) by the corresponding quantity for 
the entire catchment A C IA. This is a reasonably good ap- 
proximation when (1) A* /A* and A c /A are much less than 1 
and (2) surface runoff generation in natural catchments takes 
place predominantly by the saturation excess mechanism. 
Based on Freeze" s [1974] and Dunne's [1978] reviews of a 
large number of field studies and computer simulations, it is 
apparent that these conditions are met in most catchments, 
especially in humid areas. With this approximation, equation 
(10) simplifies to 


A * 

nXr) = Q s (t) + 0,{1 — Q s (t)] (II) 

A c 


where Q s {t) is the proportion of the total runoff rate from the 
entire catchment that is generated by the saturation excess 
mechanism and is defined by 

( 12 ) 

Note here that the approximate equation (equation (11)) is 
exact for the two extreme cases, Q s = 0 (infiltration excess 
runoff only) and Q s = 1 (saturation excess only). Also, when 
Q s = 0, I1;(t) = 0, which then yields the original GUH of 
Rodriguez-Iturbe and Valdes [1979]. 

To evaluate Ityr), we require expressions for A* /A c . In 
section 2.4 we derive expressions for A* /A c in terms of A c /A 
using geomorphologic principles and the topography-soil 
index framework of Beven [1986] and Sivapalan et ai [1987]. 


A r 

Qs(r) = ~p 

A 


A c / A c \ 

7 /> + ( 1 '7 - P' , ”' WI 


2.4. Derivation of A*JA C 

To disaggregate A c into A* for / = 1, • * • , H, we assume 
that the catchment area can be broken up into a number of 
rectangular plane hillslopes which drain directly to streams 
of different orders /. As outlined in the earlier description of 
the runoff generation model, the contributing areas are 
predicted by the topography-soil index In (aTJT x tan £). 
Given the threshold value of this index at saturation, In 
(aTJT x tan p) s , the proportion of contributing area A c /A can 
be obtained from the cumulative distribution function of In 
(aT e /T x tan 0). 

The A c /A versus In (aTJT x tan /3) relationship for an 
idealized rectangular plane hillslope was derived by Beven 
and Wood [1983]; applying this relationship to the hillslopes 
that drain directly to streams of order i and whose combined 
area as defined above is A*, we have 



where s g . is the mean ground surface slope of the hillslopes 
in a direction normal to the streams to which they drain and 
L g is the mean length of the hillslopes in the same direction. 
Horton [1945] and Morisawa [1962] have found that L g in 
most catchments can be approximated by \!2D , where D is 
the drainage density which does not vary greatly with the 
order of the catchment and can therefore be assumed to be 
constant. At very large catchment scales where there is 
significant nonstationarity in landforms, this assumption 
would not hold. 

We now make use of Horton's [1945] law of stream slopes 
which can be expressed as 

Sj/Sj - \ = Rs (14) 

where is the mean stream slope of ith-order streams. 
Strahler [1950] observed a consistent relationship between 
the mean slope of a stream and the mean ground slope of the 
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hillslope draining to it. This relationship can be approxi- 
mated by 



(15) 


where c and d are constants. Combining (14) and (15), we 
obtain a “law of hillside slopes” which can be expressed as 

VV, = ^ sc m ^5 


Substitution of (16) in (13) yields 


A c t f **/*$<; ' 

zr 



tan P SgfisG ' 


(17a) 


li-o 

A? 


a 

tan p 


Sg, R SG ' 


(17 b) 



(18) 


Combining (7), (17), and (18) and eliminating (a/tan p), we 
can derive expressions for A*IA C as a function of AJA and 
the geomorphologic constants 0, and R$g- The resulting 
expressions for a third-order catchment are presented in 
Appendix C. 

Both the variables A C IA and Q s are outputs from the 
conceptual runoff generation model described earlier. Also, 
AJA and Q, are functions of t and are dependent on the 
rainfall intensity, soil properties, topography, and the initial 
moisture conditions. The sensitivities of AJA and Q s to a 
number of dimensionless similarity parameters and auxiliary 
variables that represent climatic inputs, catchment proper- 
ties, and initial conditions were investigated by Sivapalan et 
al. [1987]. 


2.5. The Generalized GUH 

We define a dimensionless generalized GUH, denoted by 
h*(r*lr*), as 

/i*(f*lr*) = T//i(rlT) (19a) 

t» = t/Ti T*=rh l (19i>) 

where t, and t, are two characteristic time scales; r, is the 
mean duration of storms, and r, is a characteristic basin lag 
time defined here as 


TABLE I. Horton’s Order Ratios for Four Selected Catchments 


Catchment 


Santa 


North 

Bald 


Paula 

Davidson 

Nashua 

Eagle 

Parameter 

Creek 

River 

River 

Creek 

Area ratio R A 

6.8 

5.8 

5.6 

3.64 

Bifurcation ratio R B 

4.7 

4.3 

3.9 

3.17 

Length ratio R L 

1.9 

2.1 

3.2 

3.18 


Taken from Moughamian el al. [19871 and Hebson and Wood 
[1982]. 


Following Rodriguez-lturbe and Valdes [1979], /i*(/*lr*) for 
a third-order catchment is then given by 

h*(t*\r*) = B exp (- R 2 L t *) + C exp (- R L t *) 

+ (D + Et*) exp (- 2t *) (22) 

The coefficients B, C, D, and E are functions of the order 
ratios R A , R Bf R L , and R sc and the variables AJA and Q s . 
The functional relationships of 5, C, D, E to A C IA and Q s are 
expressed through the probabilities Il ( *. The expressions for 

C, D, and Eare presented in Appendix D. In the original 
formulation of Rodriguez-lturbe and Valdes [1979], n t - are 
equated to 0,-. 

The generalized dimensionless GUH was determined for 
four different combinations of the parameters R A , R B * an d 
R l and for different values of the variables A C IA and Q s . The 
values of R A , R B , and R L were taken from four catchments 
analyzed by Moughamian et al. [1987] and Hebson and 
Wood [1982]. These are summarized in Table 1. R S g was 
arbitrarily assumed to be 0.62. The variations of h p , the peak 
of the dimensionless GUH, with variations of AJA and Q s 
are presented in Figures 1 a and 1 b. They show that partial 
area runoff generation can have a significant impact on the 
GUH peak and time to peak in many catchments. 

3. Estimation of Dimensionless Flood Frequency 

3.1. Derivation of Dimensionless Peak Discharge 

Henderson [1963] has found that as long as the unit 
hydrograph peak and the time to peak are preserved, a 
triangular unit hydrograph is sufficient for the prediction of 
streamflow response. Wood and Hebson [1986] derived the 
peak discharge Q p from a catchment due to storm runoff 
generated at a constant rate q and having duration t, by 
utilizing a triangular IUH having peak h p and time to peak t p . 
The resulting expression for Q p is as follows: 

Qp = h p t r ^ 2 (23a) 


r / = Lftfv (20) 

where L ft is the length of the highest-order stream in the 
catchment and v is the mean velocity of flow in the stream 
network. Following Wood and Hebson [1986], we use an 
empirical relationship for r/ (in hours) in terms of the 
catchment area A (in square kilometers) as follows: 

r, = 2.51 A 0 38 (21) 


Q p = q Vr > 2 (23« 

In the present paper the rate of runoff generation varies in 
time even though the generated rainstorms were temporally 
constant. The results of Seven [1986] suggest that the 
temporal variations in rainfall can influence the shape of the 
flood frequency curve. This issue will be addressed in the 
discussions following the results. To capture the effects of 
the temporal variation in runoff production, the peak dis- 
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(a) 




Fig. 1. Variation of hZ, the peak of the dimensionless generalized 
GUH, to (a) A c /A and ( b ) Q s . 


charge is calculated using a dimensionless S* curve derived 
from (22). That is, 

= “7 [1-0 — exp (-/?£/*)] 
r l 

C 

+ — [1.0 - exp (~R L t*)] + 0.25(2D + E) 
r l 

• [1.0 - exp (-2r*)] - 0 .5Et* exp (-2/*) (24) 

The time-varying runoff production was then convoluted 
with the dimensionless S*(t*) curve to yield the outflow 
hydrograph from which the peak discharge and time to peak 
were extracted. As discussed in section 2, generalized GUH 
is a function of A c /A and Q s which are also time variable. 
For the simulation results presented in this paper, we used 
time-averaged values, A C !A and Q st in the calculation of h*. 

Relationships between the scaled (or dimensionless) and 
nonscaled discharge values can be established using the time 
scales T/ and r r defined earlier and the following transfor- 
mations. These transformations are similar to the ones used 
by Sivapalan et al. [1987] and are consistent with the 
definitions presented in Appendix A. 


Q P r, 


<Pc(0, - 6 r ) 

(25) 

<?T r 


- 6r) 

(26) 


The above expressions for the dimensionless peak discharge 
do not include a base flow component; its treatment is 
discussed in section 3.3. 

3.2. Frequency Distributions of Inputs 

As in most previous work, we assume that the scaled 
storm duration is exponentially distributed; in this case it has 
mean 1 . The scaled point rainfall intensities are assumed to 
be gamma distributed with parameters a and 0 and coeffi- 
cient of variation C v . Following Wood and Hebson [1986], 
we assume that the scaled mean rainfall p* over a catchment 
of area A is also gamma distributed with parameters a A and 
P A . Wood and Hebson [1986] have shown that 

a A = aK~ 2 

Pa = P* 2 (27) 

where k 2 is a geoclimatic scaling parameter defined by 

<t\ = <t 2 k 2 (28) 

where <j 2 and cr\ are the variances of point rainfall and areal 
average rainfall, respectively. The k 2 can be estimated for 
any catchment area A using the space correlation of the 
rainfall intensity process. 

The variance of point rainfall intensities within a catch- 
ment is also important and has a significant effect on predic- 
tions of the rates of runoff generation. The model of Siva- 
palan et al. [1987] requires the coefficient of variation of 
point rainfall intensities within the catchment. The relation- 
ship between and C„ for a catchment of area A is given 
by 

£^ = ( 1 _ k 2) , /2 (29) 

Following Sivapalan and Wood [1987], the spatial corre- 
lation structure of rainfall is assumed to be of the form 

Pp(r ) = a x exp (~b 2 r 2 ) + a 2 exp (-b 2 r 2 ) (30) 

where a { , a 2 , b j , and b 2 are constants with a j + a 2 = 1 . The 
correlation length is given by 



In this paper we assume a, = 0.7, a 2 = 0.3, and bfb 2 = 3.0. 
For such a correlation structure the variation of k 2 and 
C vp lC v with the scaled catchment area A* = A/A 2 is pre- 
sented in Figure 2. 

3.3. Distribution of the Initial Condition 

In the model of Sivapalan et al. [1987] the initial moisture 
state of the catchment is uniquely represented by the auxil- 
iary variable Q *. Q* is a dimensionless baseflow parameter 
which in combination with the topography-soil index param- 
eter <t> * and the hydrogeologic parameter determines the 

initial contributing areas and the soil moisture in the unsat- 
urated zone (note that the baseflow is inversely proportional 
to Q*). 

In reality, Q* can vary between storms (short term) as well 
as between seasons (long term). Ideally, the distribution of 
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Q* can only be obtained by the modeling of the processes in 
the interstorm periods. Although this is the long-term objec- 
tive of much current hydrologic research, in this paper we 
arbitrarily assume that Q* belongs to a gamma distribution 
whose parameters (mean Q* and coefficient of variation C v q ) 
can be obtained from the analysis of a large number of base 
flow recession curves of the catchment of interest. 

The base flow component of peak discharge can be related 
to the antecedent moisture conditions which then yields 




k q 


(1 + cir) w /i 


exp 



(32) 


where /, =/exp (A) and/ 2 = In the present work we 
have set/, = 400.0 and/ 2 = 0.40. Here, /is the parameter 
that controls the decrease in hydraulic conductivity with 
depth. 


scaling parameter k 2 , both of which are important parame- 
ters of the flood response. 

Realizations of the random variates are then input into the 
storm response model which results in the output Q* p . When 
arranged in ascending order, this yields the sample cumula- 
tive distribution of flood peaks P(q* p ). For n storm events per 
year the annual exceedance series flood frequency return 
periods are obtained from 

UQf) = MnP{Q*) Q* s /"'(l In ) (34) 

In the simulations carried out in this work was have used 
n = 20 storms per year and simulated a total of 4000 storms; 
this set was repeated 25 times for a total of 4000(25) storm 
analyses. For many climates, 20 storms per year may be too 
small. For the object of exploring the effect of catchment 
characterizations on the shape of the flood frequency curve, 
20 is probably sufficient. In section 4 of this paper, further 
comments are provided on how the results from this analysis 
can provide insights into the flood frequency characteriza- 
tion in actual catchments. 

4. Results and Discussion 

A number of experiments were performed with the model 
to study the sensitivity of the flood frequency response to 
many of the dimensionless similarity parameters. Some of 
these parameters (namely, p*, ^c* 0*' an ^ see 

Appendix A for their definitions) arose from the runoff 
generation model and were identified by Sivapalan et al. 
[1987]. The remainder, A*, w* 9 C^, and the Horton order 
ratios (/? x , R B > R L * and R SG )> were introduced in this paper 
during the development of the flood frequency model. The 
parameters were varied, one at a time, from the following 
arbitrary base values: p* = 0.657, Kq = 0.657, = 0.05, Q* 


3.4. Dimensionless Flood Frequency Distribution 

The cumulative distribution function for flood peaks Q* is 
given by 

F(Qf) = j f(p *, if, Q*) dp* dtfdQ* (33) 

where /(p*. f*, Q*) is the joint probability density function 
of the point rainfall intensity p*. rainfall duration, t* 9 and 
base flow parameter Q*, respectively, and R is the region 
containing all combinations of p*, /*, and Q* for which flood 
peaks exceed Q p . It is not feasible to analytically evaluate 
(33) without recourse to further simplifications. For the 
purposes of this paper, we decided to numerically evaluate 
the integral in (33) by Monte Carlo simulation. The variables 
which are assumed to vary from storm to storm are the 
scaled mean storm intensity p *, scaled storm duration t*, 
and the scaled initial soil dryness represented by Q* . As 
discussed below, the distributional characteristics of these 
variables are assumed fixed for each analysis, thus repre- 
senting a stable and stationary climate. Nevertheless, due to 
the nonlinearities in the storm response, both in the trans- 
formation from rainfall to runoff generation and in the 
relationship between flood peak and basin lag, the derived 
flood frequency distributions are not scale-independent. The 
catchment area influences basin lag t ; and geoclimatic 



Fig. 3. Variability of the flood frequency distribution over 25 
realizations. 
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Fig. 4. Value of selected parameters of the runoff generation model corresponding to the flood discharges of Figure 
3, plotted against the discharge return period. The parameters were averaged across the 25 realizations at each return 
period: (a) scaled storm intensity divided by scaled soil conductivity p*lR J, ( b ) scaled storm duration /*, (c) initial soil 
wetness Q*, ( d ) average contributing area, and (e) fraction of runoff as saturation excess runoff. 


= 0.555, <f>* = 2.34, A* = 9.0, <u* T = 1.45, R A = 6.8, R B = 
4.7, R l = 1.9, R sg = 0.62. The spatial coefficient of 
variation of the rainstorms, C v , was set at a value of 2.0 for 
all storms; for the soil parameter Kq , its C vK > was set equal 
to 1 .0, and for the initial catchment dryness its C v q was set 
equal to 1.0. The Brooks-Corey parameter B was set to 0.40, 
/i = 400.0 and / 2 = 0.40 in all simulations. Results using 
these parameter values will be referred to as the base case. 


Figures 3 and 4, presents some results for the base case. In 
Figure 3 the variability of the derived flood frequency curve 
for the 25 repetitions is shown. The average flood frequency 
curve is darkened and is essentially a straight line showing 
EV-1 behavior; the same as the rainfall-soil input. This can 
be seen in Figure 4 a where the average frequency curve for 
P*IRq is presented. To construct this, and all remaining 
figures in the paper, the 25 repetitions were averaged for 
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each return period. In Figures 4b- 4 e the frequency curves 
for scaled storm duration /*, initial catchment wetness Q*, 
average contributing area during the storm, and fraction of 
runoff due to saturated excess runoff mechanism are pre- 
sented. 

These five curves represent the significant hydrologic 
processes represented in the flood frequency curve. What is 



Fig. 5. The sensitivity of the flood frequency distribution to scaled 
storm intensity p*. 


extremely surprising is that it appears that the shape of the 
flood frequency curve is determined from the shape of the 
rainfall-soil distribution and that variables such as average 
contributing area, initial catchment wetness, and storm 
duration are constant across return periods. Further sensi- 
tivity analyses carried out later in this paper show a rather 
more complex situation than that represented here in the 
base case. Closer inspection of Figure 4 shows that the 
floods of importance to the base case catchment are heavily 
influenced by infiltration excess runoff production (Figure 
4^); so perhaps the base case represents a semiarid climate. 
In this case it is quite reasonable that initial catchment 
. wetness is independent of return period. Somewhat more 
surprising is to see that on average storm duration is only 
slightly related to return period and that contributing area is 
extremely stable at just under 30%. The constant contribut- 
ing area and increase in infiltration excess runoff with return 
period demonstrates that for the base case, catchment to- 
pography plays only a minor role in the generation of floods 
of interest. 

A series of sensitivity analyses were performed by varying 
one parameter at a time from the base case. The purpose of 
these sensitivities is to explore how fundamental hydrologic 
processes can influence flood frequency characteristics. 
Curves of the type presented in Figure 4 are available for 
each parameter set; essentially, only the flood frequency 
curves (averaged over the 25 repetitions) will be presented 
for each sensitivity analysis with selected results of other 
parameters to bring forth important features. 

Mean scaled storm intensity. Two additional runs were 
made to explore increases and decreases of the mean scaled 
storm intensity. In the base case the mean scaled storm 
intensity was set equal to the mean scaled soil hydraulic 



Fig. 6. Values of the scaled storm intensity divided by scaled 
soil conductivity corresponding to the flood discharges for the 
sensitivity run of p* = 0.986, plotted against the discharge return 
period (averaged over 25 realizations). 
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Fig. 7. Value of selected parameters of the runoff generation model corresponding to the flood discharges of the 

sensitivity run of p* = 0.329, plotted against the discharge return period. The parameters were averaged across the 25 

realizations at each return period: (a) scaled storm intensity divided by scaled soil conductivity p*IK* Q , ( b ) fraction of 
runoff as saturation excess runoff, (c) scaled storm duration /*, (d) average contributing area, and ( e ) initial catchment 
wetness Q*. 

conductivity. The two additional runs had values 50% larger clearly demonstrates the importance of both soil and rainfall 
and smaller than the base case. Figure 5 shows the results, characteristics in determining the basic shape of flood fre- 

For the case of p * larger than the base case the average flood quency curves. This is even more evident in the analysis of 

curve is still linear with the log of the return period but at a the second sensitivity curve in Figure 5. 
steeper slope. This change with respect to the base case is The second sensitivity curve had p* half of Kj a °d at first 
also exactly mimicked in the frequency curve of the rainfall- inspection the shape is somewhat confusing in that the curve 

soil conductivity ratio, as can be seen in Figure 6. This appears to go towards an EV-3 shape. Upon closer inspec- 
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Return Period (years) 

Fig. 7. (continued) 

tion, this parameter set provides a catchment where different 
processes dominate over various parts of the flood frequency 
curve. For return periods less than about 20 years the flood 
peaks are completely dominated by saturation excess runoff 
with a virtual absence of any infiltration excess runoff. At 
high return periods (greater than 100 years) the flood fre- 
quency curve reflects the frequency curve of the rainfall-soil 



Fig. 8. The envelope, over 25 repetitions, of the flood frequency 
distribution for the sensitivity run, p* = 0.329. 


distribution (as shown in Figure la). The shift in the mech- 
anisms can be clearly seen in Figure lb which shows the 
fraction of generated runoff due to saturation excess. For 
low return periods, only saturation excess runoff is gener- 
ated. For frequencies between 20 and 100 years there is a 
transition between the two curves. Since the ‘‘saturation 
excess" curve (low return periods) is controlled by catch- 
ment topography, there is no real reason why its slope 
should be the same as the "infiltration excess" flood fre- 
quency curve. 

Figures 7c-7e show the frequency curves for scaled storm 
duration, average contributing area, and initial catchment 
wetness. From these results it appears that the transition 





Fig. 9. Flood frequency distributions for three watersheds [after 
Moughamian et al. % [1987]: (a) North Nashua River, ( b ) Davidson 
River, and (c) Santa Paula Creek. 
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Fig. 10. The sensitivity of the flood frequency distribution to 
scaled soil conductivity Kq. 


from the saturation excess flood curve to the infiltration 
excess curve is accompanied by storms that are of long 
duration on a wet catchment. The uncertainty of the occur- 
rence of this combination in actual catchments contributes to 
the uncertainty as to when this transition would take place. 



Figure 8 shows the envelope for the 25 repetitions using p* = 
0.329. The transition appears to start at a scaled peak 
discharge of about 1.5, with a stable (topographically con- 
trolled) saturation excess curve and a somewhat unstable 
transition and infiltration excess flood curve. The envelope 
of Figure 8 covers a significant portion of observed flood 
data and may help explain why flexible flood frequency 
curves such as the Wakeby fit data so well. 

For actual vegetated catchments, with medium to deep 
soils having high conductivities (due to forest litter and 
macropores) when compared to mean storm intensities, such 
transitions can have an important effect on the flood fre- 
quency curve and affect the parameterization of standard 
flood frequency curves. Figure 9 presents the plotting posi- 
tions for flood data on three catchments presented by 
Moughamian et al. [1987]; it is left to the reader to see 
(imagine?) the transition described here. It should be noted 
here that the location of the transition can vary as is seen in 
the sensitivity runs for the mean scaled soil conductivity 
which are presented next. 

Mean scaled soil conductivity. In a manner similar to the 
analyses with mean scaled storm intensity, the mean scaled 
soil hydraulic conductivity was varied 50% larger and 50% 
smaller than the base case. Figure 10 presents the three flood 
frequency curves (again averaged over 25 repetitions.) The 
results essentially parallel the results for the storm intensity 
analysis. Two points should be noticed: similar changes in 
conductivity produce less sensitivity in the scaled peak flood 
discharge and the transition from saturation excess to infil- 
tration excess occurred at lower return periods for the 
parameter sets used. In Figure le there was a hint that as the 
shift to infiltration excess dominated flood curve occurred, 
the influence of the catchment soil wetness decreased (be- 



Return Period (years) Return Period (years) 

Fig. II. Value of selected parameters of the runoff generation model corresponding to the flood discharges of the 
sensitivity run of * 0.986, plotted against the discharge return period. The parameters were averaged across the 25 
realizations at each return period: (a) initial soil wetness Q* and (6) fraction of runoff as saturation excess runoff. 
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came more stable.) Figure 11a, which gives the initial 
catchment wetness for K$ = 0.986, tends to support this. 
Notice the transition at about the 2- to 5-year return period. 
Figure Mb gives the fraction of runoff due to the saturation 
excess mechanism for the same case. Notice that this 



Fig. 13. The sensitivity of the flood frequency distribution to the 
shape parameter <t>* of the distribution of In {aT e lT x tan p). 


fraction varies from 1 .0 at low return periods to about 0.40 at 
the 200 year return period. 

Initial catchment wettness. The third random input in 
the Monte Carlo simulations was the initial catchment wet- 
ness Q*. The earlier results showed that for the base case 
this parameter was stable over all return periods. Therefore 
one would not expect to see any sensitivity of the flood 
frequency curve to this parameter. Figure 12 shows this. 

Topographic and soil parameters. The flood frequency 
curve for the base case was infiltration excess dominated. 
The topographic and soil parameters, 4> * and are impor- 

tant in the saturation excess runoff generation. The hydro- 
geologic parameter ip* is a measure of the depth of the soil 
horizon and thus the depth to the water table. Small \p* 
represents a deep soil horizon and a deep water table. The 
parameter <p* is the fitting parameter for the distribution of In 
(aTJT x tan 0), the soil-topographic index introduced in 
section 2.2. Sensitivities from the base case should, there- 
fore, have little influence on the flood frequency curve. 
Figures 13 and 14 show this. Further work is being done to 
explore the sensitivity of “saturation excess” dominated 
catchments to these parameters. 

Scaled catchment area. Figure 2 presented the relation- 
ship between point rainfall variances within a catchment of 
area A for a rainfall field having a correlation length A p . The 
ratio Alkj was defined as A*> the scaled catchment area. 
Large A* correspond to “large” catchments in so much that 
the catchment is not well covered by a homogeneous rainfall 
field and the opposite for small values for A *. In the base 
case, A * was set to 9. Two additional sensitivity runs were 
performed with A* set to 1 and 36. Figure 15 presents these 
results, and it is clear that catchment scale, with respect to 
the scale of a storm, has a significant impact on the shape of 
the resulting flood frequency curve. These results support 
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Fig. 14. The sensitivity of the flood frequency distribution to the 
dimensionless soil parameter &*. 
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the results presented by Hebson and Wood [1986], which 
showed that the interplay between catchment and rainfall 
scales was critical in the shape of the derived flood fre- 
quency curve. In actual catchments the value of A* would 
vary as storms of different types occurred. It is clear from 
these results that the interplay between scales requires 
further investigation. 

5. Conclusions 

This paper has described a derived flood frequency model 
using a generalized GUH and based on partial area runoff 
generation. The first objective of the work was to obtain a 
greater understanding of the interrelationships among pro- 
cesses that underlie the storm response of catchments of 
different scales and physical characteristics as reflected in 
flood frequency distributions. We attempted to do this here 
by casting the storm response model using scaled (or dimen- 
sionless) parameters and focusing on concepts of hydrologic 
similarity. Second, we wished to understand why previous 
derived flood frequency models based on the GUH have 
performed poorly and to make significant fundamental im- 
provements to these methods. 

In this paper, we adopted a derived flood frequency 
approach with a GUH-based runoff routing model, following 
on previous work by Hebson and Wood [1982], Diaz - 
Granados et at. [1984], and Wood and Hebson [1986]. 
Significant improvements were made to these methods 
through the following extensions: (1) use of a physically 
based runoff generation model that incorporates runoff gen- 
eration on partial areas by both infiltration excess and 
saturation excess mechanisms, (2) use of a generalized GUH 
based on partial area generation and consistent with exten- 
sion 1 above, (3) variability of antecedent moisture condi- 
tions between storms, and (4) incorporation of the effects of 
catchment scale both on the rainfall input distributions and 
in runoff generation. 

It has been known for at least 50 years that the important 
problem in surface hydrology is determining "what to 
route” not "how to route." Most of the hydrology literature 
has focused on the latter topic. Our work demonstrates the 
importance of the former and the mechanism that generated 
the direct runoff. For example, the results of the Monte 
Carlo simulations have shown that for catchments domi- 
nated by infiltration excess runoff the flood frequency curve 
is completely defined by the distribution of the scaled 
rainfall-soil parameter p*!K%, and the scaled catchment area 
A*. This emphasizes the need for further research into 
rainfall distributions, especially the distributions resulting 
from storms of different types and scales. 

For catchments where saturation excess storm production 
dominates at low flood return periods and infiltration excess 
dominates at high return periods, the results show that the 
resulting flood frequency distribution may appear toward an 
extreme value type 3 (EV-3) curve, implying a limiting flood. 
Such an interpretation is incorrect, and the flood frequency 
curve is transitioning to an infiltration excess dominated 
flood curve. Further research is required to understand the 
rainfall and catchment characteristics that define the extent 
of the two mechanisms. The results of the simulation imply 
that the transition part of the flood frequency curve is 
produced by long storms of medium intensity on initially wet 
catchments. 



Fig. 15. The sensitivity of the flood frequency distribution to the 
scaled catchment area A*. 

The analysis of catchment similarity by means of the kinds 
of sensitivity analyses carried out in this work have yielded 
valuable insights into the interrelationships between various 
processes at the catchment scale. These can help in the 
development of simple physically based models of catch- 
ment behavior. The work in this paper can be extended in a 
number of directions. The first extension would be to aug- 
ment the rainfall model to include temporal variability as 
well as the spatial variability that was considered. Further, 
the Monte Carlo simulations should randomize storm scale 
(through A*) and correlate storm scale with storm charac- 
teristics to consider cyclonic and convective storms. The 
second extension is to expand the range of the Monte Carlo 
simulations by starting the base case where the mixed 
mechanisms occurred (Figure 8). This would allow one to 
explore the sensitivities with respect to those soil and 
catchment parameters which were not interesting for the 
base case presented here (e.g., Figures 12-14). 

Finally, the flood frequency model developed in this paper 
needs to be applied to some actual catchments before one 
can be sure that the conclusions made in this paper correctly 
explain the flood data presented in Figure 9. 

Appendix A 

Topography . It is assumed that In (aTJT x tan p) follows 
a three-parameter gamma distribution with location param- 
eter pi., scale parameter **, and shape parameter <f>* and 
having mean A = p + <t>*x *• 

Soil. The soil hydraulic properties are expressed in 
terms of the Brooks-Corey relationship with parameter d r , 
A,, and B [see Sivapalan et a/., [1987]. 

Hydrogeology. Q(0) is the base flow from the catchment 
at the beginning of the storm and Q 0 and / are catchment 
hydrogeologic parameters and can be estimated by the 
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analysis of the baseflow recession curve just prior to the 
storm. 

The following dimensionless parameters can then be de- 
fined: 


P* ■ 




Q*=- 


PTr 

<M*i - Or) 
*0 T r 

+c(0s ~ 0 r ) 
4>cf 


A — /x 

in imm 

A - /x 


<t> * 

c vp , C vK , and B complete the list of dimensionless similarity 
parameters. 


Appendix B 

The variables 0, denote the proportion of the hillslopes 
that drain directly into streams of order i\ The expressions 
for 0, in terms of Horton order ratios are 

e i = R 2 B /R\ (Bl) 

Rb Rb(RU^b~2) 

02 = ; (B2) 

R a R.\(2R b - 1) 


A , 


The ratios of contributing area for different orders are 



A c (A C !A) 


Kia'a -am)' 


2 OiRsG* 


i - I 


i- 1.2,3 


Appendix D 

The relationships between the parameters for the third 
order generalized GUH (equation (22)) and n,* # the number 
of particles of runoff generated on all hillslopes draining into 


streams of order / are as follows: 

B=-Y\ x A\Rl (D1) 

c= -um 2 + n 2 A 5 )flt (D2) 

D = Il,(A 3 /?£ - 2 A a ) + n 2 (A 6 *£ - 2A 7 ) (D3) 

£= -2(n,i4j + n 2 A 6 + n 3 A g )R2 (D4) 


where 


MP^Rl- 1) 

(1 - R l ){ 2 - R 2 l ) 2 


(D5a) 


UPnRL 

(\-R L )(2-R L ) 2 


(D5fc) 


R b Ra(Rl ~ 3 Rb + 2) 
R~a~ R\(2R b -X) 


(B3) 
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Appendix C 
For the following A C IA condition 

Ac 1 ^ 

— < 0 2 + 03 “ ” — 2 OftsG 
A H SG j » 2 

The ratios of contributing area for different orders are 

A * A* A* 

n C 2 J ^ 

Af Af Aq 

For the following A C IA condition 

1 3 A 3 

02 + Oi-— 2 1 efi'sh' 
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The ratios of contributing area for different orders are 
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For the following A C IA condition 
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